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We study a functional equation whose unknown maps a euclidean 
space into the space of probabihty distributions on [0, 1]. We prove 
existence and uniqueness of its solution under suitable regularity and 
boundary conditions and we characterize solutions that are diffuse on 
[0, 1]. A canonical solution is obtained by means of a Randomly Re- 
inforced Urn with different reinforcement distributions having equal 
means. 



1. Introduction. The present work treats a particular functional equa- 
tion whose unknown maps a euclidean space into the space ^([0, 1]) of prob- 
ability distributions on [0, 1]. 

Consider two probability distributions /i and v on the interval [0,/3], with 
/3 > 0, and assume that /i and v have the same mean and this is greater 
than 0. Then, for {x, y) ranging over the subspace S = [0, oo) x [0, oo) \ (0, 0) 
of M^, define the following equation with parameters fi and u: 

(1.1) X / {G{x,y)-g{x+k,y))ii{dk)+y {g{x,y)-g{x,y+k))v{dk) = 
Jo Jo 

the unknown is the function 

g ■.§^v{[o,i]). 

Without additional constraints or requirements, equation (1.1) in its com- 
plete generality admits infinitely many solutions. For instance, any constant 
function g satisfies (1.1). We will show that under suitable regularity and 
boundary conditions, the problem described by (1.1) is well-posed in the 
sense that its solution exists, it is unique and depends continuously on the 
boundary datum. Moreover, we will also prove that the solution depends 
continuously on the parameters /x and i/ and we will characterize a class 
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of solutions Q mapping the interior of S into the subspace of probabiUty 
distributions diffuse on [0, 1]. 

A particular case of the problem considered in this paper has been studied 
in [1] where it is proved that, when the two parameters /i and u are equal, 
there exists one and only one continuous solution to (1.1) that maps the 
X-axis and y-axis borders of S in the point mass at 1 or at 0, respectively, 
and that approaches the point mass at x/{x + y) as x + y tends to infinity. 
We here extend this result to the case of different parameters and with 
the same mean, and to more general boundary conditions. These will be 
described by means of a continuous function : [0, 1] "^([0, 1]) that will 
represent the boundary datum of the problem. 

In the next section we set notation and terminology, we formally describe 
the functional equation problem and we state three results concerning its 
solution; they will be proved in the rest of the paper. Section 3 deals with 
the construction of the canonical solution to (1.1) for the special case when 
the boundary datum ip{t) is the point mass at t, for all t G [0, 1]; indeed, the 
canonical solution is obtained by means of a Randomly Reinforced Urn (see 
[1, 2, 10, 11] and references therein), which represents the original stimulus 
to study the functional equation considered in the paper. Canonical solu- 
tions are the building block for proving existence, uniqueness and regularity 
properties of the solution to the functional equation problem with a gen- 
eral boundary datum; this will be shown in Section 4. Section 5 describes 
functional equation problems whose solution maps the interior of S into the 
subspace of P([0, 1]) consisting of probability distributions with no point 
masses. The final Section 6 presents some examples. Auxiliary technical re- 
sults have been postponed to the Appendix. 

2. Problem and main results. In this section we set notation and 
terminology and we describe the functional equation problem in detail. We 
also state three main results concerning its solution; they will be proved in 
the rest of the paper. 

2.1. The Wasserstein metric for spaces of probability distributions. For 
any f3 £ (0, oo), we endow the set P([0, f3]) of probability distributions on the 
real interval [0, 13] with the 1-Wasserstein metric dw which metrizes weak 
convergence. Recall that, for ^1,^2 G '^([O)/^])) 



where and are the cumulative distribution function and the quantile 
function of ^ G ^([0, /?]), respectively (see [6] for more details). Moreover, 
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by Kantorovich-Rubinstein Theorem, 



(2.1) 



mf{E{\Xi - X2I : ~ Ci, X2 ~ 6)} 



where the infimum is taken over all joint distributions for the vector of 
random variables {Xi,X2) with marginal distributions equal to .^1 and .^2) 
respectively. The metric space {V{[0, /3]),dyy) is complete and compact. 

2.2. The setT* of parameters. For < mo < /? < 00, endow the cartesian 
product P([0,/3]) X P([0,/3]) with the Manhattan-distance 



and consider the subset 7' of couples (//, v) of probability distributions with 
support in [0, /9] having means that are equal and that are both greater than 
or equal to mo, i.e. such that 



The elements of "P will act as the parameters for the functional equation 
(1.1); note that P is a closed subset of the metric space ^([0, /?]) x ^([0, f3\) 
and therefore it is compact. 

2.3. T/ie sei C([0, 1], "^([0, 1])) of boundary data. A boundary datum 93 
is defined as a continuous map from [0, 1] to ^([0, 1]). We endow the set of 
boundary data C([0, 1],7'([0, 1])) with the sup-distance 



then (C([0, 1],'P([0, l])),doo) is a complete metric space. 

From now on, 5 will indicate the element of C([0, 1],'P([0, 1])) defined by 
setting 5{t) = 6t for t G [0, 1], where 6t denotes the point mass at t. 

2.4. The set C(S, ■p([0, 1])) where solutions are to be found. Let S = 
[0,00) X [0,00) \ (0,0) and C(S,'P([0, 1])) be the set of the continuous maps 



For n = 1, 2, ... let §„ = Sn{(2;,y) & S : x + y > l/n} and consider the 
distance between elements Gi,G2 £ C'(S„,P([0, 1])) defined by 



dM((Ml,l^l),(M2,l'2)) = dwifJ-i, fJ'2) + dw{l^l,l^2) 




doo(¥'i,</^2) = sup dw{^i{t),ip2{t)); 

iG[0,l] 



dn{Qi,Q2) = sup dw{Gi{x,y),G2{x,y)). 
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We then define a new distance d by setting, for all Gi,G2 G C{S,V{[0, 1])), 
^ir r\ 1 c^n(^i |§„, ^2 |s„) 

where Q |§„ indicates the restriction to S„ of a € C(S, ^([0, 1])). 

The distance d metrizes the uniform weak convergence in any closed sub- 
set of §1J{(0)0)} which does not contain the origin. Note that convergence 
with respect to d is equivalent to convergence with respect to all d„ of the 
corresponding restrictions. The set (C(S, 'P([0, 1])), d) is a complete met- 
ric space; we will look for elements of this space that are solutions of the 
functional equation (1.1). 

2.5. The functional equation problem. The Problem object of this paper 
is now easily stated: 
given 

G V and the boundary datum (p e C([0, l],V{[0, 1])), 

find 

(2.2a) g e CiS,Vi[0,l])) 

such that, for all {x, y) € S, 



(2.2b) 



rl3 

X I {G{x, y) - Q{x + k, y))ii{dk) 

+ y / {g{x,y)-g{x,y + k))iy{dk) 



(2.2c) g{0,y)=ipiO), 
(2.2d) gix,0) =ip{l), 

(X 
Six,y),^{ 



X + y '' x+y^oo 

2.6. Main results. Our first result states that Problem (2.2) is well-posed 
in the sense of Hadamard. 

Theorem 2.1. A solution to Problem (2.2) exists, it is unique, and it 
depends continuously on the boundary datum. 

In the rest of the paper, we denote with the unique solution to 

Problem (2.2). Theorem 2.1 will be proved first in the special case when the 
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boundary datum is the map 6. Indeed Q^^ ^-j is a canonical solution for the 
problem since, for any other boundary datum ip £ C([0, 1], P([0, 1])), we will 
show that 

(2-3) ^U = ^^(^U)' 

where 

^^■.C{S,V{[0,l]))^C{S,V{[0, 1])) 
is the linear map defined by setting, for all Q £ C{S,V{[0, 1])), 

^^ig){x,y) = f\{t)g{x,y){dt) 
Jo 

with (x,y) ranging over S. 

The second theorem concerns the continuity of the solution to Problem 
(2.2) when the parameters of the equation are let to change. Indicate with G'^ 
the set of solutions to Problem (2.2) obtained by holding fixed the boundary 
datum f G C([0, 1], ^([0, 1])) and by letting the parameters (fi, u) range over 
V: 

Theorem 2.2. For any given boundary datum ip € C([0, 1], P([0, 1])), 
the map 

from 'P to G'^ , is uniformly continuous and G'^ is compact. 

To prove Theorem 2.2 we will first show that it holds for canonical solu- 
tions, i.e. for G^, and then we will prove that the map ^'(^ is continuous. 

The third result regards a different regularity property of the solution 
to Problem (2.2), which depends on the boundary datum (p but not on 
the parameters (/i,i^). Indeed we characterize solutions 0"^^^-^ mapping the 
interior of S into the class of probability distributions on [0, 1] having no 
point masses; such solutions will be called diffuse. 

A boundary datum E (^([0, 1],7'([0, 1])) is said to be monotonic if, for 
ah s,t£ [0, 1], s<t, 

•^{s) <st V>{t)- 

For a given ip G C([0, 1], V{[0, 1])), indicate with $ the probability distribu- 
tion on [0, 1] obtained as the convex combination with uniform weights of 
the members of the family {</?(*) : t E [0, 1]}; i.e. $ = Jq ip{t)dt. 
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Theorem 2.3. Assume that the boundary datum ip is monotonic and let 



1. If there is (a^Ojyo) the interior of S such that ^.^(a^Oi yo) has no 
point masses in [0, 1], then $ = Jq (p{t)dt has no point masses in [0, 1]. 

2. If ^ = Jq ip{t)dt has no point masses in [0,1], then G^^y^{x,y) has no 
point masses in [0, 1] for all (x, y) in the interior of S. 

Once again, in Section 5, we will first prove Theorem 2.3 for canonical 
solutions and then for the general solution Q"^^ . 

3. Existence of canonical solutions: a Randomly Reinforced Urn. 

In this section we assume that the boundary datum of Problem (2.2) is the 
map 6 € C([0, 1], 7'([0, 1])); we prove the existence of a solution for this spe- 
cial instance of the problem, by constructing it through a generalized urn 
scheme. This solution will be called canonical since the solution to Problem 
(2.2) for a general boundary datum will be obtained by transforming the 
canonical solution through a suitable map. The generalized urn scheme that 
we are going to consider is the Randomly Reinforced Urn (RRU) introduced 
in [4, 5, 8, 10, 11] and further studied in [1-3, 9, 12]. This is an urn containing 
balls of different colors, say black and white. The urn is sequentially sampled; 
whenever a black ball is extracted from the urn, it is replaced in it together 
with a random number of black balls having distribution ^, whereas if the 
color of the extracted ball is white, the ball is replaced in the urn together 
with a random number of white balls having distribution u. The extra balls 
added every time the urn is sampled are called reinforcements. 

On a rich enough probability space (17, A, V), define two independent infi- 
nite sequences of random elements, {Un} and {{Vn, Wn)}; {Un} is a sequence 
of i.i.d. random variables uniformly distributed on [0,1], while {(V"„,VFn)} 
is a sequence of i.i.d. bivariate random vectors with components uniformly 
distributed on [0,1]. Then, define an infinite sequence {{Rx {n) , Ry (n))} of 
bivariate random vectors by setting, for all n. 



where and are the quantile functions of two distributions /x and i' 
having support in [0,/3], with /3 > 0. Let x and y be two non-negative real 
numbers such that x -|- y > 0. Set Xq = x, Yq = y, and, for n = 0, 1, 2, let 



be the unique solution to Problem (2.2). Then: 



Rx{n) = q^,{Vn) and Ryin) = g^(W„) 



(3.1) 



( 




Xn + Rx{n+l)l{n + l), 
y„ + i?y(n + l)(l-I(n+l)) 
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where the variable I(n + 1) is the indicator of the event {Un+i < + 
Yn)~^}. The law of {{Xn,Yn)} is that of the stochastic process counting, 
along the sampling sequence, the number of black and white balls present in 
a RRU with initial composition {x,y) and reinforcement distributions equal 
to fi and ly, respectively. 

For n = 0, 1, 2, . . . let = Xn + y„ be the total number of balls present 
in the urn at time n and set 

Zn{x,y) = Xn/Dn] 

Zn{x, y) represents the proportion of black balls in a RRU with initial com- 
position before the (n + l)-th ball is sampled from it. In [11] it is 
proved that {Zn{x,y)} is eventually a bounded sub- or super-martingale, 
and it thus converges almost surely, and in L^, for 1 < p < oo, to a random 
variable Z^{x^y) G [0,1]; moreover, when /u and u have different means, 
Zoo{x-, y) is the point mass concentrated in 1 or 0, according to whether the 
mean of ^ is greater or smaller than that of v. However, when the means of 
/i and u are the same, the distribution of Z^{x^y) is unknown, apart from 
a few special cases, see [1] and [9]. 
For a given couple (/i, v) € V, let 

/:(^,,) : S ^ P([0, 1]) 

be the map which assigns to every (x, y) G § the distribution of the limit 
proportion Zoo(x,y) of a RRU with initial composition {x,y) and reinforce- 
ment distributions /i and u. In the special case where = i^, the map jC(^^^) 
has been characterized in [1] as the unique solution to Problem (2.2) when 
the boundary datum is 5. We now extend this result to the general case 

Proposition 3.1. ^{^i^v) is a solution to Problem (2.2) when its bound- 
ary datum is equal to 5. 

In order to prove Proposition 3.1 we need some auxiliary results; when 
they do not depend on the parameters (//, v) G T', and there is no place for 
misunderstanding, we write C for jC(^,i/) . Some technicalities connected with 
the Doob's decomposition of the process {Zn} have been postponed to the 
Appendix. 

The distance between C, evaluated at {x,y), and the boundary datum, 
evaluated at x/{x -\-y), is controlled in the following lemma; this distance is 
uniformly bounded, provided that the size of the urn initial composition is 
sufficiently large. 
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Lemma 3.1. Ifx + y>2(j, 



x+y' ^ X + y 

Proof. Note that, by (2.1), dw{C{x,y),6_^J = E{\Z^{x,y) - ^|). 
Set 

m = kfi{dk) = / kv{dk)] 



then m > niQ. From Lemma A. 2 and Lemma A. 3 with n = 0, and because 



X + y > 2P imphes \/ < 1, we have 



x + y x + y yx + y yx + y 

□ 

The Markov inequahty together with Lemma 3.1 imply the following 
corollary. 

Corollary 3.1. If x + y >2j3, 

x I ^ \ ^ 2 / ^ 



x + y^ J hoV x + y 



p(|Zoo(x,y)-^— I >ho) < 
for every /iq > 0. 

Lemma 3.2. For all hq > 1 and e > 0, there is N = N{e,nQ) such that, 
E {\Zn{x,y) - Zoo(x,y)|) < e, 
if n > N and x + y > 1 /hq . 
Proof. Equation (A. 4) yields 

(3.2) P(I)„ <t)=F(^>'-)<tE±<t^^ - 



Dn t^ Dn mo{n - 1) + P 



for all t > 0. Set 

(3.3) t = max{16/3/e^2/3} 
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and 

. ^ -{l + noiP -mo)) - (3 

(3.4) N > ^ — — — - + 1. 

mo 

From (3.2) and (3.4), we get 

(3.5) FiDN <t)<l- 

Moreover, since the process {{Xn, Yn)} is Markov, it follows from Lemma 3.1 
and (3.3) that, for n > iV and lo £ {Dn > t}, 
(3.6) 

E(\Z^-Zn\ {Xn,Y^)){uj) = dw( C{Xn{uj),Y,{uj)), 6 ) 



e 

< -. 
2 



Since {-Dn+i < ^ i^n < t} for all n, (3.5) and (3.6) imply that 



E[|Zoo -Zn\]= E[|Zoo -Zn\; {Dn > t}] + E[\Z^ -Zn\;{Dn< t}] 

< E(E(|Zoo - Z„|l{B„>i}|(X„,y„))) +P(Z)jv < t) 

< e 

for n>N. □ 

The next result can be read as a bound on the modulus of continuity of 
C when evaluated at the inner points of S. 

Lemma 3.3. For all hq > 1 and e > 0, there is rj = T/(e,no), increasing 
with € and I/uq, such that 

dw{C-{x,y),C{x,y)) < e, 

if \x — x\ + \y — y\ < rj and min{x + y,x + y} > I/uq. 

Proof. Let N = A^(e/4,no) be given by Lemma 3.2. Then: 

dw{C{x,y),C{x,y)) 

< E[\Zoo{x,y) - ZN{x,y)\] + E[\Zoo{x,y) - ZN{x,y)\] 
+E[\ZN{x,y) - ZN{x,y)\] 

<i + E[\ZN{x,y)-ZN{x,y)\]- 

For controlling the last term, we adopt a coupling argument as in [1]. 
Consider two different randomly reinforced urns, the first one with initial 
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composition {x,y) and second one with {x,y). The two urns are coupled 
in the sense that the same processes {Un} and {(Kul^n)} generate both 
{{Xn{x,y),Yn{x,y))} and y), y)) according to the dynamics 

described in (3.1). With the same arguments as in [1, pages 701-702], one 
may show that 

E[\Zn{x, y) - Zn{x, m < (1 + N)- 1^ " ^1 + I?/ " ?/l - 



min{x + y, X + y} ' 
therefore, if ?? < 2{i+N)no ^ 

n\ZN{x,y)-ZM{xM] < 



□ 



Proof of Proposition 3.1. By considering the conditional distribu- 
tion of Zoo{x,y), given I(l),i?x(l) and i2y(l), and taking the expected 
values, one immediately verifies that -C^^ satisfies equation (2.2b) for all 
(x, y) € S. Conditions (2.2c) and (2. 2d) are also easily verified when ip = 5. 
Finally, (2.2a) and (2.2e) are consequences of Lemma 3.3 and of Lemma 3.1, 
respectively. □ 

The next result proves a further regularity property of ^(/j.^^) ■ 

Proposition 3.2. The map 

from {T',d][i) to (C(S, 'P([0, 1])), d), is uniformly continuous. 

Proof. Let A : C(S,P([0, 1])) xV^ C(S, P([0, 1])) be the operator 
defined by setting, for every H G C(S,P([0, 1])) and {fi,u) G V, 

A{n, (/i, i^))ix, y) = -^yj^ nix + k, y)fi{dk) + ^yj^ nx, y + k)i^{dk) 
■ n{x + q^{t), y)dt + ^J^ n{x, y + q,{t))dt, 



X 

x+y 



where (x, y) ranges over S. 
Let n> 1. Then 



(3.7) d„(A(Wi,(^,Z.))|§„,A(7^2,(/i,I^))k) <f^n(Wlk,W2k), 
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for every 7ii,7i2 G C(S,'P([0, 1])) and (fJ-,!^) € T*. Indeed, for every (x,y) G 



Moreover, if 7^ G C(S,'P([0, 1])) is Lipschitz on S.„ with Lipschitz constant 
Kn, then, for every (//i, z/i), (^2, 1^2) G "P, 

(3.8) (in(^(7t:, {fli,Ui))\s^,A{n, (/U2,Z^2))|S„) < A'nC?A/((W,'^l), (/^2,l^2)), 

since, for every (x,?/) G Sn, 

dw{A{n, {fii,ui)){x,y),A{n, {fi2,i^2)){x,y)) 



Now, for every H G C(S,P([0, 1])) and G set A^iH, = H 

and, for N = 1,2, ... define iteratively 



Consider Hq G C(S, V{[0, 1])) defined by setting 7^0(2;, y) = '5(^) for every 
(x, y) G S; then Zq^''^\x, y) has distribution Ho{x, y), while, for = 1, 2, 
Zj^''^\x,y) has distribution {Hq, {fi,u)){x,y), where, for clarity of ex- 
position, the exponent of the Z variables is evidence for the reinforcement 
distributions of the RRU under consideration. Note that, for n > 1, TCq is 
a Lipschitz map from S„ to ^([0, 1]) with Lipschitz constant n. Moreover, 
it is not difficult to show, with computations analogous to those appear- 
ing in (3.9), that the operator A preserves the Lipschitz property with the 
same constant; hence, for (/x,i/) G "P and N = 1,2,..., {TCo,{ijl,i')) is a 
Lipschitz map from S„ to ^([0, 1]) with Lipschitz constant n. 




< sup dwi'Hi{x',y'),n2ix',y')). 



(3.9) 




x + y Jo 
< K„dA/((;Ui,i/i), (^2,1^2)) • 



{n,{fi,u)) = A{A^-\n,{fi,u)),{fi,u)). 
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Let (/^i, J^i), (/X2, J^2) G "P, re, > 1 and, for ease of notation, set Tii = 
A^-\no,{^ii,Ui)), for i = 1,2; then 
(3.10) 

dniA^in^, (W, l/l ))!§„, (7^0, (/U2, 1^2)) Is J 

= dn{A{ni, (;Ui,l/i))|§„, A(H2, (/U2,1^2))|§„) 

< Ci„ (^(Hl , , I'l ) ) !§„ , ^(^2 , (m , l/l ) Is J ) 

+ (i„ (^(^2 , (m , 1^1 ) ) l§„ , ^( ^2 , (^^2 , 1^2 ) ) Is J 

< C^n(^l|§„,'H2|§„) + ndA/((m> '^l), 

the last inequality being a consequence of (3.7) and (3.8). By iteratively 
applying (3.10), it follows that 

dn{A^{no, (^^1, i^i)), A^iHo, {fi2, 1^2))) < nNdM{{fii,iyi), (/i2, 1^2)). 

Therefore, for every n > 1 and e > 0, if iV = A^(e, n) is chosen according to 
Lemma 3.2, we obtain 

dni^ilMiMi) kn,'^{f,2M2) ISn) < (^(/^i ) \s„ , A^ {Hq , {fll , 

+ dniA'^iHo, (/XI, ZVl))|§„,^^(Wo, (/U2, i^2))\n„ 
+ C?n(-4^(Wo, (/i2,'^2))|§„,'C(^2_,^2) |s„) 

< nNdM{ifJ-i,i^i), {fJ-2,'^2)) 

+ sup \E{\Zt^'^^\x,y)-Zi;'''''\x,y)\) 
(x,j/)e§„ 

+ E(|z};^'^^)(x,y)-Z(^^^'-^)(x,y)|) 

< nNdM{{p.i,yi), ifJ-2, 1^2)) + 2e. 

This shows that the map (/i, u) ^ ^{n,v) l§„5 from (7-*, ^m) to (C(S„, ^([0, 1])), (i„) 
is continuous for every n. Hence the map (/i, z^) 1— > from {V,dM) to 

(C(§, "PdO, 1])), d) is continuous; since P is compact, it is also uniformly 
continuous. □ 

4. The solution to the functional equation problem. In this sec- 
tion we prove Theorem 2.1 and Theorem 2.2. In particular we show existence 
and uniqueness of the solution to Problem (2.2) when the boundary datum 
is a generic element of C([0, 1],'P([0, 1])). Existence is shown by means of a 
constructive proof based on the canonical solution described in Section 3. 
Uniqueness is proved through a fixed point argument. 
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Given ip £ C([0, 1]V{[0, 1])), define the map : P([0, 1]) ^ V{[0, 1]) by 
setting, for every ^ € ^([0, 1]), 

rAO{B)= j\{t){B)adt), 
Jo 

where B ranges over the Borel sets in [0, 1]. 

Lemma 4.1. For any given (p € C([0, 1]7^([0, 1])), the map is uni- 
formly continuous. 

Proof. Since (p G C([0, 1], ^([0, 1])), (p is uniformly continuous and bounded: 
i.e. for any e > 0, there is an = 7/(e, p) such that 

(4.1) dw{^{ti),p{t2)) <e, if |ti -t2| < ??, 
while 

(4.2) dw{^{ti),p?{t2)) < 1, for all ti,t2 G [0,1]. 

Now, take ^i, ^2 ^ ^([0) 1]) such that dw{S,i,(,2) < ^V- We are going to prove 

that diy(r^(ei),r^(e2)) < 2e. 

Because of (2.1), there is a probability space {Q,A,W) carrying a couple 
of random variables Xi,X2 such that Xi ~ ^1, X2 ~ ^2, and dwiS,ii(,2) = 
E(|Xi — X2I) < er]] by Markov inequality, 

(4.3) F{\X^ - X2\ > 7]) < e. 

On the product probability space {0. x [0, 1], ^ (g) B{[0, 1]),F ® Leb) define 
the random variables 

Ci(w,t) = inf \^z: J^^ ^ p{r]i{u;)){ds) > t} = q^(r,i{uj)){t) 

and 

C2(w,t) = inf {z: J^^ ^ 99(r?2(u;))(c?s) > = g<^(^2(<^))(t), 

where indicates the quantile function relative to the probability distribu- 
tion ^ G V{[0, 1]). 

For i = 1,2, note that <p{Xi) is the conditional distribution of d, given 
Xi] thus, Ci ~ r<^(^j). Moreover, for all w G 

diy((/9(Xi(u;)),(^(X2(u;))) = |g^(Xi(a;))(i) - 9<^(X2H)(*)|c^* 

|Ci(w,t)-C2(c^,t)|dt. 
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Hence, 

dH'(r^(ei),r^(6)) <iE(|Ci-C2|) 

(4.4) =E(E(|Ci-C2||Xi,X2)) 

= ¥.{dw{^{Xi),^{X2))). 

Now, let F = {\Xi - X2I > r/}. From (4.1), (4.2) and (4.3) one obtains: 

E {dw[^{Xi), ^{X2))) = E {dw{^[Xi), (^(^2)); F) + E {dw{'f{Xi), (^(^2)); F'' 

< F{F) + e]P(F^) < 2e. 

The last inequality, together with (4.4), concludes the proof. □ 



Proof of Theorem 2.1.. (i) Existence. When the boundary datum ip 



6, the existence of a solution G^tn,^) is guaranteed by Proposition 3.1, and 



this is ^(^,^) ■ 

Now let if e C([0, l],V{[0, 1])) and define, for ah {x,y) € S, 

we are going to show that is indeed a solution to Problem (2.2) when 

the bondary datum is (p. In other words, the composition 

{x,y) ^ Ql^^^){x,y) ^ TviQf^,u)ix,y)) 

gives a solution to Problem (2.2), i.e., (2.3) holds if the map is defined 
by setting, for all G G C(S,P([0, 1])), 

^^iG){x,y)=T^iGix,y)) = f\{t) G{x,y){dt) 

Jo 

with {x,y) ranging over S. 

Because of Proposition 3.1, Gf^ satisfies (2.2b) when the border datum 
is 5. Since is linear, this implies that, for all (x, y) G §, 



r, 



(r^ / C^)(^ + ^'^H^^) + 7T^ / Gl,^){x,y + kMdk)) 



'^'x + yj ''^'"'"'^ - ' x + y 

^^^{Gl,.)){x + k,y)f,idk) + -^ [ ^^{Gl^,^){x,y + kMdk); 



x + y J ^^''"^'^ " ^ ' x + y 



hence ^(^(^(^j^)) satisfies (2.2b) when the border datum is ip. Now, by 
Lemma 4.1, ^(/^(^(^iy)) is a continuous map from S to ^([0,1]), being the 
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composition of the continuous maps and Q^^j^yy hence (2.2a), (2.2c), 
(2. 2d) and (2.2e) also hold true. 

(ii) Uniqueness. Sketch of the argument. Our argument in [1, Section 5] 
can be easily extended to this more general situation. 

Condition (2.2e) requires Q to be continuous at the projective infinite 
points. It is therefore convenient to transform the space S along the projec- 
tive automorphism r of so defined: 

{x : y : u) ^{u : x : x + y). 

The automorphism r has the following properties: 

- the space S is mapped into the affine space S* = [0, oo) x [0, 1]; 

- the positive x-axis is mapped into itself by (x,0) — > (l/x,0); 

- the positive y-axis is mapped into the semifine {y* = l,x* > 0} by 
(0,y) ^ (l/y,l); 

- the projective infinite point relative to the direction = kis mapped 
in the point (0, k); 

- the origin is mapped in the projective infinite point (1:0: 0). 

The inverse map of r is (x* : y* : u*) [y* : u* — y* : x*). Problem (2.2) 
can be equivalently formulated on S* as follows: 

given 

{lJL,v) G V and the boundary datum e C([0, 1],7'([0, 1])), 

find 

(4.5a) g* ^C{S\V{[Q,l])) 

such that, for all (x,y) € S*, 



(4.5b) 



(4.5c) r(x*,0) =(/j(0), 

(4.5d) g*{x\l)=^{l), 
(4.5e) g*{Q,y*) = ^{y*). 
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In fact, (4.5b) is just (2.2b) in the new coordinates. Indeed, the transforma- 
tions 



Gix,y) = g*{T{x,y)) 



g*{x\y*) 



if (x*,y*) G (0,(») X [0,1]; 



lim e?(r-i(s*,y*)) if x* = 0, y* G [0, 1], 



show the equivalence of Problem (2.2) and Problem (4.5). 

Now, let C* (S*) be the space of continuous function Ti* : S* ^ V{[0, 1]) 
such that, for every {x*,y*) G S*, 

n*{x*, 0) = ifiO), n*{x\ 1) = ^{l) and 7^*(0, y*) = ip{y*). 

Define the following operator A* mapping C* (S*) into C* (S*): 

u{dk) 



A*{n*){x*,y*) = y* J n*[j 



+ ii-y*) J n*[- 



+ kx* ' 1 + kx* 
x* y 



+ kx* 1 + kx 



with (x*,y*) ranging over S*. 

With the same argument used in [1, Theorem 5.2], one can prove that A* 
has at most one fixed point; hence Problem (4.5) has at most one solution. 

(iii) Continuity with respect to the boundary datum. We prove this last 
part by showing that 

(4-6) d(gf;^),gf;,))<doo(^i,^2) 

for all ipi,ip2 G C([0, 1],P([0, 1])). We recall here that (see, e.g., [6]) 



(4.7) dw{vi,m) = sup{| J h{t)r]i{dt) - J h{t)r]2{dt) 

where is the Lipschitz norm. Then, for h such that 

{x, y) G S, we get 



L < 1 and 



/ h{s)gl;,^^{x,y){ds) 



< 



h{s) J ^i{t){ds)gl^,){x,y){dt) 

-J h{s) J Mt)ids)gl,,)ix,y){dt) 
h{s)Mt)ids) - [ h{s)Mt)ids) gl,,){x,y){dt) 



< / dw{^i{t),(p2{t))g(^^^i,){x,y){dt) < doo{^i,<f2); 
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hence dw{Q'(^„){x,y),g'^^^^^{x,y)) < (ioo(v?i, again by (4.7). Inequality 
(4.6) follows easily. □ 

Remark 4.1. Given € T', the inequality (4.6) can he completed 

as follows: for all <^i,(/72 £ ^^([0, 1],'P([0, 1]), 

Hence, for any given {li^v) G "P, we have an embedding 

C([0,1],P([0,1]))^C(S,P([0,1])). 
Indeed, for n = 1, 2, (2.2e) implies that 

since dn < 1, and thus dn < 2 -^^"^ , ^/lis implies the right inequality in (4.8). 

Remark 4.2. Let m he the common mean of (/-f,i^) G "P. For p € 
[mo/m, 1], sei (/x',i^') = (pyt^ + (1 -15)501?"^ + (1 -p)'^o)- Then {i-i',^') G P 
and 

Remark 4.3. Let h : [0, 1] — > [0, 1] be a continuous function and ip G 
C([0, 1], P([0, 1])) a boundary datum. For ^ G 7^([0, 1]), denote with /lo^ the 
distribution of the random variable h{W), where W is a random variable 
with probability distribution ^. Then 

for all (x, y) G S. 

Remark 4.4. One may notice that the boundary conditions (2.2c) and 
(2. 2d) are redundant. Indeed, if (2.2b) and (2.2e) are true for a Q : S ^ 
7^([0, 1]), then Q satisfies (2.2c) and (2. 2d). To see this, let x > and 
consider Q{x,0). By iteratively applying (2.2b), one obtains 



I 

13 rP 



JO 

P rP 



Q{x + ki+ k2, 0)iJ.{dki)fi{dk2) 

/ •••/ g{x + ki^ \- kn,0)fi{dki) ■ ■ ■ fj.{dkn) 

Jo Jo 



G. ALETTI,C. MAY,P. SECCHI/A FUNCTIONAL EQUATION P{[0, 1]) VALUED18 

for all n = 1,2, .... However, because of (2.2e), ifJ2i'=i h ^ oo as n —>■ oo, 
then 

n 

lim dw{Q{x + y ki, 0), V3(l)) = 0. 

i=\ 

Hence, the Law of Large Numbers and the Dominated Convergence Theorem 
imply that 

g{x,o) = ^{i). 

The argument for proving that Q{0,y) = (p{0), if y > 0, is analogous. 

We are finally in the position to prove Theorem 2.2. 

Proof of Theorem 2.2.. The theorem is true when the boundary da- 
tum is 5, as follows from Proposition 3.2 and the fact that T' is compact. 
For a general boundary datum ip G C([0, 1],P([0, 1])) the result follows once 
it is proved that the map is continuous; but this is a consequence of 
Lemma 4.1. □ 

5. Diffuse solutions. As an immediate consequence of Proposition 3.1 
and [2, Theorem 3.2] we have the following result, which is a particular 
instance of Theorem 2.3 and represents a tool for proving it. 

Proposition 5.1. For all {x,y) in the interior ofS and z € [0,1], 

Sl,u)i^,y)i{z})=o. 

Proof of Theorem 2.3. Given z € [0,1], note that (fz{t) = (f{t){{z}) 
is a measurable function of t, since it is the monotone limit of the sequence 
of continuous functions fci"^ defined by setting, for n = 1, 2, ... and t G [0, 1], 

4")(t) = I hi-HsMt){ds), 

where 

The two functions 

ip~{t) = sup (w), '^t{t) = i^<^(t) {z) 



ns — nz + 1, \{ z — \/n < s < z, 
nz — us -\- 1, \i z < s < z + 1/n, 
0, otherwise. 
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are monotonically nonincreasing in t, since is monotone. Moreover, (pz{t) = 
— ip~{t). Therefore is a bounded variation function and it thus has 
at most a countable number of points of discontinuity. Note that 

^({4)= t Ut)dt= [\ci^t{t)-<t^-At))dt- 

Jo Jo 



Proof of part 1. Let (a;o,yo) be a point in the interior of S such that ^^(xq, yo) 
has no point masses in [0,1]. By way of contradiction, suppose there is a 
z £ [0, 1] such that ^>({2}) > 0. 

Since <l>({z}) > 0, there are e > 0, a > and G [0, 1] such that 

(5.1) ^z{t) > e, for all t € /* = [z* - o, + a] n [0, 1]. 

Set 

64/?, 



R= {(x,y) e S: X > max(2/?, ■ 



a2 



y=(=-l)x (z=— ^ 
z X + y 

a a 

For all {x,y) € R, Corollary 3.1 with /lo = f implies that 



ze[z^- -,z* + ^]ri [o,i]|. 



p(Zoo(x,y) 0/,) <P(jZoo(x,y)-z| > ^) < -J^ < I 
V /V 2/ava;2 

and thus Of^ ^){x,y){h) > |. Then, because of (5.1), for all {x,y) € R, 
Sl,.)i^,y){{^}) = f\.{t)Gl^,){x,y){dt) 



Now, set 

r = inf{n > 0: {Xn{xo,yo),Yn{xo,yo)) £ R} (inf{0} = oo); 

it is not difficult to show that P(r < oo) = p > 0. Thus the strong Markov 
property implies that Gff,^^){xo,yo){I*) > |; therefore ^^^^^-,(xo, yo)({^}) > 
p| contradicting the assumption that G'^^ i,^{xo,yo) has no point masses in 
[0, 1]. This concludes the proof of part 1. 
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Proof of part 2. Let now z € [0, 1]; by assumption $({z}) = 0. Since ipz has at 
most a countable number of points of discontinuity, (/^^ > and / ipz{t)dt = 
0, there is a sequence ti,t2i ••• such that ipz{t) = for all t € F = (Ujjtj})'^. 
Then, given any (x, y) in the interior of S, 



J [0,1] 

vz{t) Ol,u)i^, y){dt) + ^ Mti) y){{ti}) 

[o,i]nF ^ 
= 

the last term being zero because of Proposition 5.1. □ 

6. Examples. In this section we give explicit descriptions of the so- 
lution for some specific choices of the reinforcement distributions 
(/i,z>) G V and of the boundary datum if G C([0, 1], 7^([0, 1])). The first 
example is prototypical since it considers the Polya urn scheme and the 
family of Beta distributions, whose properties had a central role in originat- 
ing most of the problems tackled in this paper. 

6.1. The Polya urn scheme and the family of Beta distributions. We 
indicate with Beta(2;, y) the beta distribution on [0, 1] with parameters 
(x, y) £ S. If (x, y) is a point in the interior of S, Beta(x, y) has a den- 
sity given by 

f r(x -I- y) x-i/-, +\y-i 

/Beta{x,j/)(ij = Y{x)r{y) ^ ~ ' ' 

for t G [0, 1]; it is convenient to indicate with Beta(0,y) and Beta(x,0) the 
point mass at or at 1, respectively. 

The random limit composition of a Polya urn with initial composition 
(x,y) G S and constant reinforcement equal to 1 has distribution Beta(x,y); 
indeed a Polya urn is a special RRU with reinforcements ^ = u = 6i. The 
same result holds when /x = z/ = Bernoulli(p), for p > 0; see [4] and the 
references therein. Hence, for p > 0, 

= ^{Bcrnoulli(p),Bernoum(p))(^'y) = Beta(x,y), 

for all (x, y) G S. 

The Polya urn scheme, where both reinforcement distributions are equal 
to the same point mass and the boundary datum is 5, may be easily extended 
by considering a more general boundary datum along the hint given in Re- 
mark 4.3. For instance, let 7 > 0, define y?^ : [0,1] — > 7'([0, 1]) by setting 
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f'yit) = 6{t^^"'), for t G [0, 1], and consider Problem (2.2) with fi = ly = 5i 
and boundary datum equal to 99^. Note that (p^ is monotone and thus the 
unique solution to the problem is diffuse. Indeed for (x, y) in the interior S, 
the distribution Q^^^ g^^{x,y) has density 

^{Sj,6-^)y-^^y> T{x)T{y) 

for t G [0, 1]. For x = 1 and y > 0, the solution QJ^^ ^^^{l,y) is called the 
Kumaraswami distribution with shape parameters 7 and y; see [7]. 

6.2. Bernoulli reinforcements. A more intriguing extension of the Polya 
urn scheme is obtained by considering reinforcement distributions {fJ^,i^) G 7' 
different from equal point masses; we here treat the case where /i and u are 
scaled Bernoulli distributions with the same mean. Let m > kfj^ > > 0, 
and assume that and v are the distributions of two random variables, 
say Rx and Ry, such that Rx/kfj, has distribution Bernoulli(m/A;^) while 
Ry /kn has distribution Bernoulli(m/A;j,). 

Equation (1.1), with (^,1^) as above, reads 

which does not depend on m. One easily verifies that the equation is satisfied 
by the continuous map Q : S ^ "^([0, 1]) defined by setting, for all (x, y) G S, 

X y 

Q{x,y) = Beta(— , — ). 

Moreover, note that, 

rf^(Beta(f ,f),^( Z^^- )) 0. 
fc^ ky xky + yk^ {x+y)^oo 

Hence, if h : [0, 1] — > [0, 1] is defined by setting 

tki, 



h{t) 



tky + {l-t) 



for all t G [0, 1], then 



gf;J)(x,y) = Beta(- , , 
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for {x,y) € S, is the unique solution to Problem (2.2) when fi and v are the 
scaled Bernoulli distributions defined above and the boundary datum is the 
continuous map h o 5 : [0, 1] — > ^([0, 1]) defined by setting 

h o 6it) = 5{ ^^"^ 



tk^ + {l-t)k^' 
for all t £ [0,1]. 

We now want to find the distribution of the limit composition of a RR U 
whose reinforcements are distributed according to the scaled Bernoulli dis- 
tributions fi and z^. Note that h is continuous, monotonically increasing and 
its inverse is 

UKfj_ + (1 — U)Ki, 

for u € [0, 1]. Then it follows from Remark 4.3 that, 

Gl,.){x, y) = G\-^l:{'°'{x, y) = h'' o Beta(^, |-), 

for all (x, y) G S. For (x, y) in the interior of §, u)i^' v) ^ density and 
this is 

tol,^^i.,y)it) - k. k^ r(x/Mr(2/A.)[,,^ + (i_,)fe^]t+^ 

for t £ [0,1]. 

APPENDIX A: DOOB DECOMPOSITION OF THE RRU PROCESS 

This Appendix provides a series of auxiliary results necessary to prove 
Propositions 3.1 and 3.2. We will refer to the notations introduced in Section 
3. For n = 1, 2, ... let An = cr{6i, Rx{l), Ry{^), ■ ■ ■ , ^n-, Rxin) , Ry {n)) and 
consider the filtration {^n}; then, given the initial urn composition (x,y) € 
S, the Doob's semi-martingale decomposition of Z.„(x,y) is 

Zn{x, y) = Zo{x, y) + M„(x, y) + A„(x, y) 

where {M„} is a zero mean martingale and the previsible process {An\ is 
eventually increasing (decreasing), again by [11, Theorem 2]. We also denote 
by {(M)„} the bracket process associated to {M„}, i.e. the previsible process 
obtained by the Doob's decomposition of M^. 
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We first provide some auxiliary inequalities. As a consequence of [2, 
Lemma 4.1], we can bound the increments A An of the Z„-compensator pro- 
cess and the increments A(M)„ of the bracket process associated to {M„}. 
In fact, an easy computation gives 



and 
where 

and 



A^+l = E{AZn+l\An) = Zn{l " ^n)^;+l 



E{{AZn+lf\An) = Zn{l " ^n)^:+l. 



•^n ) ) 



Rx(n+l) Ryin+l) 

A* Tff I Dn Dn 

- L , Rx{n+l) , Rvin+l) 

Dn Dn 



Rx{n+l) 2 RY{n+l) 2 



Zl^.=^({l-Zn){ i^Vj ^H±Tl) K 

Now, [2, Lemma 4.2] with m = Jq kfi^dk) = Jq ki'{dk) gives 

By applying [2, Lemma 4.1] with h{x,t) = (^)^ Bd = [2/3, oo), D = Dn, 
R = Rx{n + I) ov R = RY{n + 1) and A = An-, one obtains: 



(A.2) Z„(l-Z„)(-^j <E((AZ„+i)^|A)<^n(l-^n;^^^^^^2' 

on the set {L>„ > 2/3}. Since 

E((AZ„+i) Vn) = E((AA+i + A/W„+i) Vn) = {AAn+if + A(M)„+i, 
if Dq > 2/3, and thus (3 + Dn> 3/3, (A.l) together with (A.2) yields 

m \2/ / (3 — m 



m/3 



A{M)n+l < Zn{l ~ Zn 



{P + DnY' 
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Lemma A. 1. For all k = 1,2, 

1 + (/3 - m)/Do 



(A.4) 



1 



< 



'Dk' - Do + m{k-l)+P 
If, in addition, Dq > 2/3 then, for all k,n = 1, 2, 



(A.5) 



E 



1 



1 



when d > c> and bk = c + Dn — /3 + rnk. 



< 



(3-m + d / I 



1 



m 



Proof. Let r]* be a random variabile independent of and let rji 
be a random variable independent of cr{Aoo,V*) and such that has 
distribution Binomial(l, Define -4.^_,_„- = (j{ri* , Ak+n-i,%k + n)); by 

[2, Lemma 4.1], if D > is -measurable and < R< P with E{R) = 

m is independent of -4.^,_,_„- , one obtains 



E 

and thus 
(A.6) E 



1 



D + R 



A* 



1 



< 



m 1 



+ 



f3 — m 1 



PD + 13 (3 D 



E 



1 



D + 7JI 



A* 

■^k+n- 



Dk+n + V* 



A 



I(/c + n)E 



1 



Dk+n-i + V* + Rx{k + n) 



A* 



+ (l-I(A; + n))E( 



1 



Dk+n-i + if + RY{k + n) 



A 



k+n 



< I{k + n)K 



1 



^Dk+n~i +r]* + r]i 
+ (1 -I(A; + n))E 



At 

1 



E 



1 



^Dk+n-i +r]* +r]i 



Dk+n~i + if + m 
A1 



A 



k+n 



Therefore, for c > 0, by applying (A.6) A;-times, we get 



(A.7) 



E 



1 



Dk+n + C 



An] < E 



1 



Dn + c + rjk 



An 



where 7]k is independent of a{Aoo) and has distribution Binomial(A;, m//3). 
Equation (A.4) is now a consequence of [13, Eq. (21)]: if rfk ~ Binomial(A;, r) 
and I > 0, 



1 



1 
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Apply this to (A. 7) with 7i = 0, rfk = 7]^/ P, I = Dq/ /3 and r = m/ (3 to obtain 
(A.4). 

Equation (A. 5) is a consequence of [13, Eq. (25)]: if 1% ^ Bino(A;,r) and 
/ > 1, 



< 



^l + rfk'' l + kr-{l-r) 

Apply this to (A. 7) with rfk = Vk/P, I = Dn + c/f3 (which is greater than 2) 
and r = m/P to obtain 



E 



1 



An] < 



1 



c + Dn + m{k + l) - p' 



Jensen's inequality yields E,{{d + Dn+k) ^\A.n) > {d + Dn + mk) ^, and thus 



E 



1 



An 



< 



P — m + d — c 



{c + Dn + m{k + 1) - P){d + Dn + mk) ' 



Since 



m 



bk bk+i {c + Dn- P + mk){c + Dn- P + mik + 1)) 
we get (A. 5): 



E ( -rn- 



1 



+k d+Dn+k 



1 1_ 

bk bk+i 



< 



< 



P — m + d — CC + Dn — /3 + mnk 
m d + Dn + mk 

P — m + d — c 
m 



□ 



The following Lemma A. 2 and Lemma A. 3 provide inequalities which 
control the previsible and the martingale part of the process Z„ respectively; 
they require that the initial composition of the urn is sufficiently large. 



Lemma A. 2. If Dq > 2/3, then 



E(sup|A^|) < 



Dn 



Proof. Apply (A. 5) with n = 0, c = m, d = p. Equation (A.l) then 
reads 

E(|^U,|)<(2^-,„)(1--L). 



k+l' 
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if 6/; = km + Dq — {(3 — m). Since Aq = 0, 

E(sup I Al) < E ( E l^^fc+il) < E I n\Al+A) 
^ k k 

1 \ _2[5-m 



4 bk+i^ 4 Do-{P-m) 

<JL, 

where the last inequahty is true beacause (3 — m < (3 < Dq/2. □ 
Lemma A. 3. Let Dq > 2/3. For all n > 0, 

E((M)oo - {M)n\An) < 

Proof. Since Z„,+fc(l - Zn+k) < 1/4, by (A. 3), one gets 

A , , ,v m / 1 1 

A{M)n+k+i < TnTT^ vi < 



4(/3 + ) 2 - 4 V Z)„^_fc p + 

Apply (A. 5) with c = and d = (3, obtaining 

„-,/A/,rv i<N m2/? — m/1 In 
E(A(A/),,+fc+i|A) < -J— i-r- - 7— , 

if hk = km + Dn — (3. Thus 

'~ A. 



E((M)oo - {M)n\An) = E ( E A(M)fe+„+i 

k>0 

2(3 -m / 1 1 
— ~A 2^ 



4 ^^pV6^ 

< 1 < A 

- 2 2(Z)„-/3) - L»o' 

since 2(D„ - /?) > 2(Do - /3) > Dq. □ 
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